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In this work we study the dynamics of universe in f(R) = R 2 — Rq modified gravity with Palatini 
formalism. We use data from recent observations as Supernova Type la (SNIa) Gold sample and 
Supernova Legacy Survey (SNLS) data, size of baryonic acoustic peak from Sloan Digital Sky Survey 
(SDSS), the position of the acoustic peak from the CMB observations and large scale structure 
formation (LSS) from the 2dFGRS survey to put constraint on the parameters of the model. To 
check the consistency of this action, we compare the age of old cosmological objects with the age of 
universe. In the combined analysis with the all the observations, we find the parameters of model 
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as R = 6.192Io;i?7 x H$ and a 
PACS numbers: 04.50. +h, 95.36,+x, 98.80.-k 
I. INTRODUCTION 



Recent Observations of Supernova type la (SNIa) pro- 
vide the main evidence for the accelerating expansion of 
the Universe [1,2]. Analysis of SNIa and the Cosmic Mi- 
crowave Background radiation (CMB) observations in- 
dicates that about 70% of the total energy of the Uni- 
verse is made by the dark energy and the rest of it is in 
the form of dark matter with a few percent of Baryonic 
matter [3-5]. The "cosmological constant" is a possible 
explanation of the present dynamics of the universe [6]. 
This term in Einstein field equations can be regarded as 
a fluid with the equation of state of w = — 1. However, 
there are two problems with the cosmological constant, 
namely the fine-tuning and the cosmic coincidence. In 
the framework of quantum field theory, the vacuum ex- 
pectation value is 123 order of magnitude larger than the 
observed value of 10 -47 GeV 4 . The absence of a funda- 
mental mechanism which sets the cosmological constant 
to zero or to a very small value is the cosmological con- 
stant problem. The second problem known as the cosmic 
coincidence, states that why are the energy densities of 
dark energy and dark matter nearly equal today? 

There are various solutions for this problem as the de- 
caying cosmological constant models. A non-dissipativc 
minimally coupled scalar field, so-called Quintessence can 
play the role of this time varying cosmological constant 
[7-9]. The ratio of energy density of this field to the 
matter density in this model increases by the expansion 
of the universe and after a while dark energy becomes the 
dominant term of the energy-momentum tensor. One of 
the features of this model is the variation of equation 
of state during the expansion of the universe. Various 
Quintessence models like k-csscncc [10], tachyonic mat- 
ter [11], Phantom [12,13] and Chaplygin gas [14] provide 
various equations of states for the dark energy [13,15-21]. 

Another approach dealing with this problem is using 



the modified gravity by changing the Einstein-Hilbert ac- 
tion. Recently a great attention has been devoted to this 
era because of the prediction of early and late time ac- 
celerations in these models [22]. While it seems that the 
modified gravity and dark energy models are completely 
different approaches to explain cosmic acceleration, it is 
possible to unify them in one formalism [23]. 

In this work we obtain the dynamics of the modified 
gravity f(R) = \J R 2 — Rq in the Palatini formalism [24] 
and use the cosmological observations as SNIa, SDSS 
acoustic peck in CMB and structure formation to put 
constraint on the parameter of the action. In Sec. II we 
derive the dynamics of Hubble parameter and scale fac- 
tor in Palatini formalism for general case of f(R) gravity. 
In Sec. Ill using action f(R) = \/ R 2 — Rq we obtain the 
dynamics of universe. In Sec. IV we use the observa- 
tions from the evolution of background as SNIa, CMB 
and baryonic acoustic oscillation to constrain the param- 
eters of the model. Sec. V studies constraints from the 
large scale structure formation and in Sec. VI we com- 
pare the age of universe from the model with the age of 
old cosmological objects. The conclusion is presented in 
Sec. VII. 



II. MODIFIED GRAVITY MODELS IN PALATINI 

An alternative approach dealing with the accelera- 
tion problem of the Universe is changing the gravity law 
through the modification of action of gravity by means 
of using f(R) instead of the Einstein-Hilbert action. For 
an arbitrary action of the gravity there are two main ap- 
proaches to extract the field equations. The first one is 
the so-called metric formalism in which the variation of 
action is performed with respect to the metric. In the 
second approach, Palatini formalism, the connection and 
metric are considered independent of each other and we 



have to do variation for those two parameters indepen- 
dently. General form of the action in palatini formalism 
is : 

S\f;g,t,iSf m ] = ~ J d 4 x^gf(R) + S m [g M „,tf m ], 

(1) 

where k = 8ttG and S^g^u, ^/ m ] is the matter action 
which depends only on metric g^ u and on the matter 
fields * m . R = R(g,f) = g^R^(f) is the generalized 
Ricci scalar and R^ v is the Ricci tensor of the affine con- 
nection which is independent of the metric. Varying the 
action with respect to the metric results in: 



f(R)R fW (f) - l -f{R)g^ = kT„ v , 



(2) 



where prime is the differential with respect to Ricci scalar 
and T^y is the energy momentum tensor 



-2 6S m 



(3) 



Varying the action with respect to the connection and 
after contraction gives us the equation that determines 
the generalized connection as: 



V a [/'(ii)V=fl«r] = o, 



(4) 



where V is the covariant derivative with respect to the 
affine connection. We can see that the connections are 
the Christoffel symbols of the new metric h^ v where it 
is conformally related to the original one via the h^ v — 

Equation (2) shows that in contrast to the metric vari- 
ation approach, the field equations are second order in 
this formalism which seems more intuitive to have sec- 
ond order equations instead of fourth order. On the other 
hand forth order differential equations have the instabil- 
ity problem [28]. We use FRW metric for universes as 
follows 



ds 2 = -di 2 + a(t) 2 S ij dx i dx j 



(5) 



and assume universe is filled with perfect fluid with the 
energy-momentum tensor of T" = diag(-~p,p,p,p), tak- 
ing the trace of equation (2) gives, 



Rf'(R) - 2f(R) = kT, 



(0) 



where T = g^T^ = —p + 3p. Using Generalized Ein- 
stein equation (2) and FRW metric we obtain generalized 
FRW equation as [29]: 



(H !/>_ l «(p + 3p) 1 / 
2/' 6 /' 6/'' 



(7) 



equation we can express the time derivative of Ricci 
scalar as: 



R = +3H 



(l-3p')(p + p) 
Rf"-f'(R) 



(8) 



Using equation (6) we obtain the density of matter in 
terms of Ricci scalar as: 



Kp 



2/ - Rf 
1-3cj ' 



(9) 



where w = p/p. Substituting (9) in (7) we obtain the 
dynamics of universe in terms of Ricci scalar: 



H 2 = 



1 



3(l+w)/-(l + 3o;)i?/ / 



6(1 - 3w)/' 



, il-x , , ,\ f'(Rf'-2f) 



(10) 



On the other hand using equation (6) and continuity 
equation, the scale factor can be obtained in terms of 
Ricci scalar 



1 



Kpo(l — 3a>) 



(2/ - Rf) 



3(1+") 



(11) 



where po is the energy density at the present time and we 
set ao = 1. Now for a modified gravity action, omitting 
Ricci scalar in favor of the scale factor between equations 
(10) and (11) we can obtain the dynamics of universe (i.e. 
H = H(a)). For the matter dominant epoch u> = 0, these 
equations reduce to: 



1 



3/ - Rf 



6/' 



3 f"(Rf'-2f) 

2 fw-n 



— {V-Rf) 

npo 



(12) 
(13) 



III. MODIFIED GRAVITY WITH THE ACTION 

OF F(R) = sjR 2 - Rl 

In this section our aim is to apply the action of f{R) = 
\J R? — Rq in equations (10) and (11) to obtain the dy- 
namics of universe. For simplicity we use dimensionless 
parameters in our calculation defined by 



f(R) = H 2 JX*-X t 



(14) 



in which Hq is the Hubble parameter at the present time 



and X = and Xq 



II;: 



For convenience wc can 



write the action and its derivatives as: 



f{R) = H 2 F(X) 



(15) 



For the cosmic fluid with the equation of state of p = 
p(p), using the continuity equation and the trace of field 



f(R) = F'(X) 



(16) 



f"(R) 



F"{X) 



(17) 



where derivative at the right hand side of equations is 
with respect to X (i.e. ' = -rg)- We can write (12) with 
new parameter of X as: 

1 3F - XF' 

H -y X )--^- 3 F''(XF'-2F)-> 



6F' (1 



2 F>(XF"-F<) 

where H(X) 

at the present time. Using the definition of tt m (X) = 
KPm/{3H 2 ) from equation (9) we can obtain Cl m (X) in 
terms X as: 

, ( x) = ?*zX£L. (19) 



H 2 /Hq and #o is the Hubble parameter 



\1 r, 



We use fl m (X) = S} m a~ 3 (where Q m = Cl-m) an d substi- 
tute F in terms of X from (14) to obtain scale factor in 
terms of X as: 

-1/3 



a = 



1 



30, 



X 2 - 2X 2 



Vx 2 



(20) 



where to have positive scale factor, X should change in 
the range of X > V2X - It should be noted that this 
model has only one free parameter of Xq- X p represents 
the value of X at the present time and from (18) we can 
find X p in terms of Xq. On the other hand substituting 
X p = g{Xo) in (19) we will have direct relation between 
Xq and Cl m . So knowing Xq from observation one can 
calculate also Cl m . 

In the next section we will compare observations with 
the dynamics provided by this model in the matter dom- 
inant epoch. However we can see how the universe ex- 
pands at the radiation dominant epoch. We set p = l/3p 
which implies traceless energy momentum tensor and us- 
ing our proposed action in equation (6), we get a constant 
Ricci scalar of R = \/2Rq for this area. Substituting this 
constant curvature in equation (7) we have: 

(21) 

Since i?o is in the order of Hq (sec [24]), for the early 
universe we can ignore second term at the right hand side 
of this equation which results the scale factor to increase 
as a oc t 1 / 2 . 



IV. OBSERVATIONAL CONSTRAINTS FROM 
THE BACKGROUND EVOLUTION 

In this section we compare the new SNIa Gold sam- 
ple and supernova Legacy Survey data, the location of 
baryonic acoustic oscillation peak from the SDSS and the 
location of acoustic peak from the CMB observation to 
constrain the parameters of the model. We choose vari- 
ous priors applied in the this analysis as shown in Table 
I 



TABLE I. Different priors on the parameter space, used in 
the likelihood analysis. 



Parameter 


Prior 




Qtot 


1.00 


Fixed 




0.020 ± 0.005 


Top hat (BBN) [25] 


h 




Free [26,27] 


w 





Fixed 



A. Examining the Model by Supernova Type la: 
Gold Sample 

The Supernova Type la experiments provided the main 
evidence for the present acceleration of the universe. 
Since 1995 two teams of the High-Z Supernova Search 
and the Supernova Cosmology Project have discovered 
several type la supernovas at the high redshifts [17,30]. 
Riess et al.(2004) announced the discovery of 16 type la 
supernova with the Hubble Space Telescope. This new 
sample includes 6 of the 7 most distant (z > 1.25) type 
la supernovas. They determined the luminosity distance 
of these supernovas and with the previously reported al- 
gorithms, obtained a uniform 157 Gold sample of type 
la supernovas [31-33]. Recently a new data set of Gold 
sample with smaller systematic error containing 156 Su- 
pernova la has been released [34]. In this work we use 
this data set as new Gold sample SNIa. 

More recently, the SNLS collaboration released the 
first year data of its planned five-year Supernova Legacy 
Survey [35]. An important aspect to be emphasized on 
the SNLS data is that they seem to be in a better agree- 
ment with WMAP results than the Gold sample [36] . 

We calculate the apparent magnitude from the f(R) 
modified gravity and compare with new SNIa Gold sam- 
ple and SNLS data set. The supernova measured ap- 
parent magnitude m includes reddening, K correction, 
etc, which here all these effect have been removed. The 
apparent magnitude is related to the (dimensionless) lu- 
minosity distance, Di,, of an object at redshift z through: 



M + 5logD L (z;X Q ), 



where 



Also 



D L (z;X ) = (l + z) 



o H(z') 



(22) 



(23) 



M = M+5l0g {r l ik?) +25 - (24) 

where M is the absolute magnitude. The distance mod- 
ulus, ^, is defined as: 




FIG. 1. Distance modulus of the SNIa new Gold sam- 
ple in terms of redshift. Solid line shows the best fit 
values with the corresponding parameters of h = 0.63, 



0.276 



+0.376 



X = 6.207 



+0.230 
0.147 



with x 2 min/N d . .f = 0.912 for f(R) model. 



in lcr level of confidence 




FIG. 2. Distance modulus of the SNLS supernova data in 
terms of redshift. Solid line shows the best fit values with 
the corresponding parameters of h — 0.70, Q. m = 0.246_t ' 1 g 4 
and Xq = 6.484_toggg in lcr level of confidence with 
X 2 rmn/N d .a. S = 11.766 for f (R) model. 



M 



51og D L (z;X ) 

c/H 



-5 log 



1 Mpc 



25, 



[i = 5\ogD L (z:X ) + M. 



(25) 



(26) 



To compare the theoretical results with the observational 
data, we calculate theoretical distance modulus. Dis- 
tance modulus can be written in terms of new parameter 
X as: 

_1{XF'-2F)<& f x XF"-F' dX 
L ~3 (3f2 m )l Jx p (XF -2F)iH(Xy ( ' 

To put constraint on the model parameter, the first 
step is to compute the quality of the fitting through the 
least squared fitting quantity y 2 defined by: 



X 2 (M,X ) 



(28) 



where Oi is all of observational uncertainty in the distance 
modulus. To constrain the parameters of model, we use 
the Likelihood statistical analysis: 



(29) 



where Af is a normalization factor. The parameter M is 
a nuisance parameter and should be marginalized (inte- 
grated out) leading to a new \ 2 defined as: 



X = -2 In 



+ 00 



(30) 



Using equations (28) and (30), we find: 



X 2 (X ) = x 2 (M = 0,X ) - +ln(C/27r), (31) 



where 



B(X ) = 



C 



[Hobs{zi) ~ ^th{zi\ X , M = 0)] 



and 



(32) 



(33) 



Equivalent to marginalization is the minimization with 
respect to M. One can show that x 2 can be expanded in 
terms of M as [37] : 

XsNia(*o) = X 2 {M = 0, X ) - 2MB + M 2 C, (34) 

which has a minimum for M = B/C: 

B(X a ) 2 



XsNIa(^0)=X 2 (M = 0,X ) 



C 



(35) 



Using equation (35) we can find the best fit val- 
ues of model parameters, minimizing XsNia(^o)- Us- 
ing new Gold sample SNIa, the best fit values for the 
free parameter of the model is Xq = 6.207io'i4? states 



Q m = 0.276±g;|S, with Xrmn/Nd.o.f = 0.912 at la level 
of confidence. The corresponding value for the Hubble 
parameter at the minimized \ 2 is h — 0.64 and since 
we have already marginalized over this parameter we do 
not assign an error bar for it. The best fit values for 
the parameters of model by using SNLS supernova data 



is X = 6.484 



+0.099 
-0.099 



indicates Q, 



0.246 



+0.164 
-0.164 



with 



Xmin/Nd.o.f = 11-766 at la level of confidence. Figures 
1 and 2 show the comparison of the theoretical predic- 
tion of distance modulus by using the best fit value of 
X and observational values from new Gold sample and 
SNLS supernova la, respectively. In figures 5, 6, relative 
likelihood for Xq are indicated. We report the best value 
of Xo at 1 and 2-er confidence level and other derived 
parameters in table II. 



B. CMBR Shift parameter 

Before last scattering, the photons and baryons are 
tightly coupled by Compton scattering and behave as 
a fluid. The oscillations of this fluid, occurring as a 
result of the balance between the gravitational interac- 
tions and the photon pressure, lead to the familiar spec- 
trum of peaks and troughs in the averaged temperature 
anisotropy spectrum which we measure today. The odd 
peaks correspond to maximum compression of the fluid, 
the even ones to rarefaction [38] . In an idealized model of 
the fluid, there is an analytic relation for the location of 
the m-th peak: l m ks ml a [39,40] where I a is the acoustic 
scale which may be calculated analytically and depends 
on both pre- and post-recombination physics as well as 
the geometry of the universe. The acoustic scale corre- 
sponds to the Jeans length of photon-baryon structures 
at the last scattering surface some ~ 379 Kyr after the 
Big Bang [5]. The apparent angular size of acoustic peak 
can be obtained by dividing the comoving size of sound 
horizon at the decoupling epoch r s (zd ec ) by the comoving 
distance of observer to the last scattering surface r(zdec)'- 



Tj_ _ r s {z dec ) 

I A ~ r(zdec) ' 



(36) 



The size of sound horizon at the numerator of equation 
(36) corresponds to the distance that a perturbation of 
pressure can travel from the beginning of universe up to 
the last scattering surface and is given by: 



° v s (z')dz' 
ca H(z')/H 



(37) 



where v s (z)~ 2 = 3 + 9/4 x pb(z) / ' p ra d{z) is the sound 
velocity in the unit of speed of light from the big bang 
up to the last scattering surface [19,39] and the rcdshift 
of the last scattering surface, z dec , is given by [39]: 



z dec = 1048 [1 + 0.00124(u> b )-°- 738 ] [1 + <?iK>) 92 ] 



,9i = 0.0783M" - 238 [1 + 39.5(uj b ) - 763 ] 
. 92 = 0.560 [l + 21.lM 1 - 81 r\ 



(38) 



where uj m = fl m h 2 and o>b = Vl^ti 2 . Changing the param- 
eters of the model can change the size of apparent acous- 
tic peak and subsequently the position of I a = t^/Oa 
in the power spectrum of temperature fluctuations on 
CMB. The simple relation l m w ml a however does not 
hold very well for the first peak although it is better for 
higher peaks [2]. Driving effects from the decay of the 
gravitational potential as well as contributions from the 
Doppler shift of the oscillating fluid introduce a shift in 
the spectrum. A good parameterization for the location 
of the peaks and troughs is given by [40,41] 



lm =Ia(™>— <Pm) 



(39) 



where <f> m is phase shift determined predominantly by 
pre-recombination physics, and arc independent of the 
geometry of the Universe. Instead of the peak locations 
of power spectrum of CMB, one can use another model 
independent parameter which is so-called shift parameter 
1Z as: 



1Z oc 



'i 



flat 



h 



(40) 



where l{ lat corresponds to flat pure-CDM model with 
Q m = 1.0 and the same uj m and uj as the original model. 
The location of first acoustic peak can be determined in 
model by equation (39) with <f>\{u) m ,u)b) ~ 0.27 [40,41]. 
It is easily shown that shift parameter is as follows [42] : 



- Dl (Zdec'i Xo) 
(1 + z dec ) 



(41) 



In writing the shift parameter in this form we have 
implicitly assumed that photons follow geodesies deter- 
mined by the Lcvi-Civita connection. Furthermore, in 
order to use the shift parameter, the evolution of the uni- 
verse is considered in such a way that at the decoupling 
we have standard matter dominated behavior. Since we 
do not have an explicit expression for the Hubble pa- 
rameters in terms of redshift, it is useful to rewrite the 
shift parameter in terms of dimensionless parameter X 
as follows: 



K = xl£l m H 2 



H(z) 
a'(R) dR 



Rdec a 2 (R) H(R) 

F' - XF" dX 
x p (2F-XF')iH(X) 



(42) 



The observational results of CMB experiments corre- 
spond to a shift parameter of TZ = 1.716±0.062 (given by 



-i— i — |ii — |ii — i— I—i — i— r- 

Flat X = 5.0 0- m = 0.396 
Flat X = 6.0 D- m = 0.298 
Flat X = 7.0 flm = 0.187 

lcdm Q m = 0.396 , £1 = 0.604 
lcdm £l m = 0.298 , £l x = 0.702 
lcdm Cl m = 0.187 , Q. x = 0.813 
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FIG. 3. Evolution of density contrast in the f{R) gravity 
model versus redshift for different values of Xq. For compar- 
ison we plot 5 for ACDM with the same value of fi m .(A.U. 
stands for arbitery unit) 



WMAP, CBI, ACBAR) [5,43]. One of the advantages of 
using the parameter 1Z is that it is independent of Hubble 
constant. In order to put constraint on the model from 
CMB, wc compare the observed shift parameter with that 
of model using likelihood statistic as [42] : 



C 



Xcmb/^ 



where 



2 _ [7^-ob 

XCMB — ' 



n 



the 



2 

°CMB 



(43) 



(44) 



C. Baryon Acoustic Oscillations 

The large scale correlation function measured from 
46, 748 Luminous Red Galaxies (LRG) spectroscopic 
sample of the SDSS (Sloan Digital Sky Survey) includes 
a clear peak at about 100 Mpc [44]. This peak was 
identified with the expanding spherical wave of baryonic 
perturbations originating from acoustic oscillations at re- 
combination. The comoving scale of this shell at recom- 
bination is about 150Mpc in radius. A dimensionless and 
independent of Hq of this observation is: 



A 



HoD L (z a dss', Xq) 



H(z sdss ; X )z 2 dss (l + z sdss ) 2 



1/3 



(45) 
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FIG. 4. Growth index versus redshift for different values of 
Xo- To make sense we plot / for ACDM with the same Q m . 





f Zsdss dz 


_ ^sdss J 


o H(X)\ 



(46) 



we can write the above dimensionless quantity in term of 
our model parameter as: 



"(3a»n 



3z sd 



(F' - XF")dX 



H(X)(2F - XF') 2 / 3 



(47) 



We can put the robust constraint on the f(R) modified 
gravity model using the value of A = 0.469 ± 0.017 from 
the LRG observation at z s d ss = 0.35 [44]. This observa- 
tion permits the addition of one more term in the \ 2 °f 
equations (35) and (44) to be minimized with respect to 
H(z) model parameters. This term is: 



Asdss 



\A 



obs 



A 



the 



(48) 



SDSS 



This is the third observational constraint for our analysis. 



D. Combined analysis: SNIa+CMB+SDSS 

In this section we combine SNIa data (from SNIa new 
Gold sample and SNLS), CMB data from the WMAP 
with recently observed baryonic peak from the SDSS to 
constrain the parameter of modified gravity model by 
minimizing the combined x 2 = XsNia + Xcmb + xIdss 
[45]. 

The best values of the model parameters from the fit- 
ting with the corresponding error bars from the likelihood 




FIG. 5. Marginalized likelihood functions of f(R) modi- 
fied gravity model free parameter, Xq. The solid line corre- 
sponds to the likelihood function of fitting the model with 
SNIa data (new Gold sample), the dashdot line with the 
joint SNIa+CMB+SDSS data and dashed line corresponds 
to SNIa+CMB+SDSS+LSS. The intersections of the curves 
with the horizontal solid and dashed lines give the bounds 
with lcr and 2a level of confidence respectively. 




FIG. 6. Marginalized likelihood functions of f(R) mod- 
ified gravity model free parameter, Xq. The solid line 
corresponds to the likelihood function of fitting the model 
with SNIa data (SNLS), the dashdot line with the joint 
SNIa+CMB+SDSS data and dashed line corresponds to 
SNIa+CMB+SDSS+LSS. The intersections of the curves 
with the horizontal solid and dashed lines give the bounds 
with la and 2a level of confidence respectively. 



function marginalizing over the Hubble parameter in the 
multidimensional parameter space are: Xq = 6.1691q '^4 
corresponds to Q m = 0.28lloJsI at ^ a confidence level 
with Xmin/Nd.o.f = 0.902. The Hubble parameter corre- 
sponding to the minimum value of x 2 is ft, = 0.63. Here 
we obtain an age of 14.56^' 28 Gyr for the universe. Us- 
ing the SNLS data, the best fit values of model parameter 
is: X = 6.428±g; 92 states Q m = 0.2521°;^ at lcr confi- 
dence level with Xmin/Nd.o.f = 11-582. Table II indicates 
the best fit values for the cosmological parameters with 
one and two a level of confidence. Relative likelihood 
analysis for Xq using CMB and SDSS observations are 
shown in figures 5 and 6. 



V. CONSTRAINTS BY LARGE SCALE 
STRUCTURE 

So far we have only considered observational results re- 
lated to the background evolution. In this section using 
the linear approximation of structure formation we ob- 
tain the growth index of structures and compare it with 
the result of observations by the 2-degree Field Galaxy 
Rcdshift Survey (2dFGRS). The continuity and modified 
Poisson equations for the density contrast 8 = Sp/p in 
the cosmic fluid provide the evolution of density contrast 
in the linear approximation (i.e. 5 <C 1) [46,47] as: 



8 + 2-8- k 2 V 2 + 4nGp] 8 = 0, 
a L J 



(49) 



the dot denotes the derivative with respect to time. The 
effect of modified gravity in the evolution of the struc- 



tures in this equation enters through its influence on the 
expansion rate. 

The validity of this linear Newtonian approach is re- 
stricted to perturbations on the sub-horizon scales but 
large enough where structure formation is still in the lin- 
ear regime [46-48] . For the perturbations larger than the 
Jeans length, Xj = n^Vg/y/Gp, equation (49) for cold 
dark matter (CDM) reduces to: 



8 + 2-8 

a 



AnGpS = 0. 



(50) 



The equation for the evolution of density contrast can be 
rewritten in terms of the scale factor as: 



da 2 



da 



2H 

a 



3HZ 
2d 2 a 



;Sl m 5 = 0. 



(51) 



In order to use the constraint from large scale structure 
we should rewrite above equation in terms of X. So we 
have 



aH(X), 

aH 2 (X) + a 2 H(X)H'(X) 



dX 
da 



(52) 



Using equation (52), equation (51) becomes: 



dH 
da 2 



dS_ 

da 



n'(x) dx 

H(X) da 



3fi r 



2H 2 (X)a ! 



■8 = (53) 



Numerical solution of equation (53) in FRW universe 
for the two cases of f(R) gravity and ACDM model with 
the same matter density content is compared in Figure 
(3). 



In the linear perturbation theory, the peculiar velocity 
field v is determined by the density contrast [46,49] as: 



v(x)=^/*(y)-^ 



rd 3 y, 



where the growth index / is defined by: 



din 5 
din a' 



(55) 



and it is proportional to the ratio of the second term of 
equation (50) (friction) to the third (Poisson) term. 

We use the evolution of the density contrast S to com- 
pute the growth index of structure /, which is an impor- 
tant quantity for the interpretation of peculiar velocities 
of galaxies [49,50]. Replacing the density contrast with 
the growth index in equation (51) results in the evolution 
of growth index as: 



df 
din 



a 



3H O _ f2 _ f 

2a 2 a 



1 



a 



H 2 



(56) 



The above equation in terms of dimcnsionlcss quantity 
is: 



df 

din a 2aH 2 (X) 



r-f 



2 + 



aH'(X) dX 



H{X) da 



(57) 



Figure (4) shows the numerical solution of equation 
(57) in terms of redshift. As we expect increasing Xq 
causes decreasing Q m from best fit and a decreasing in 
the evolution of density contrast versus scale factor and 
growth index in the small redshifts. This behavior is the 
same as what happens in the ACDM model. 

To put constraint on the model using large structure 
formation, we rely to the observation of 220, 000 galax- 
ies with the 2dFGRS experiment provides the numerical 
value of growth index [44]. By measurements of two- 
point correlation function, the 2dFGRS team reported 
the redshift distortion parameter of /3 = f/b = 0.49±0.09 
at z = 0.15, where b is the bias parameter describing the 
difference in the distribution of galaxies and their masses. 
Verde et al. (2003) used the bispectrum of 2dFGRS 
galaxies [51,52] and obtained b ver d e = 1.04 ± 0.11 which 
gave / = 0.51 ± 0.10. Now we fit the growth index at 
the present time derived from the equation (57) with the 
observational value. This fitting gives a less constraint 
to the parameters of the model, so in order to have a 
better confinement of the parameters, we combine this 
fitting with those of SNIa+CMB+SDSS which have been 
discussed in the previous section. We perform the least 
square fitting by minimizing X 2 = x! N ia+XcMB+XsDSS+ 
Xlss. where 

2 [fobs(z = 0.15) -fth(z = 0.15; X )} 2 
Xlss — ~2 ^°°) 

The best fit value with the corresponding error bar 
for the Xq by using new Gold sample data is: Xq = 



0.278J278 73 at lcr confidence 



6.19218;$? provides fi„ 
level with xL in / N d.o.f = 0.900. Using the SNLS super- 
nova data, the best fit value for model parameter is: Xq = 

at lcr confidence level 



(54) 6.433±g;g|? gives Vt m = 0.252 



+0.147 
-0.150 



with Xmin/Nd.o.f = 11-486. The error bars have been 
obtained through the likelihood functions (£oc e~ x '/ 2 ) 
marginalized over the nuisance parameter h. The likeli- 
hood functions for the four cases of (i) fitting model with 
Supernova data, (ii) combined analysis with the two ex- 
periments of SNIa+CMB, iii) combined analysis with the 
three experiments of SNIa+CMB+SDSS and (iv) com- 
bining all four experiments of SNIa+CMB+SDSS+LSS 
are shown in Figures 5 and 6. The best fit value and 
age of universe computing in the f(R) modified gravity 
model are reported in table II. 



VI. AGE OF UNIVERSE 

The age of universe integrated from the big bang up 
to now for flat universe in terms of Xq is given by: 



to(X ) 



to 



dt 



dz 



F'-XF" dX 
3H J Xo 2F - XF> H(X) 



(59) 



Figure 7 shows the dependence of H to (Hubble param- 
eter times the age of universe) on Xq for a flat universe. 
In the lower panel we show it for ACDM for compari- 
son. As we expect modified gravity behaves as a dark 
energy and increasing Xq result in a longer age for 
the universe in the f(R) modified gravity model. 

The " age crisis" is one the main reasons of the accel- 
eration phase of the universe. The problem is that the 
universe's age in the Cold Dark Matter (CDM) universe 
is less than the age of old stars in it. Studies on the old 
stars [53] suggest an age of 131^ Gyr for the universe. 
Richer et. al. [54] and Hasen et. al. [55] also proposed 
an age of 12.7 ± 0.7 Gyr, using the white dwarf cooling 
sequence method (for full review of the cosmic age see 
[5]). To do another consistency test, we compare the age 
of universe derived from this model with the age of old 
stars and Old High Redshift Galaxies (OHRG) in var- 
ious redshifts. Table II shows that the age of universe 
from the combined analysis of SNIa+CMB+SDSS+LSS 

is 14.6918H Gvr and 13 - 61 io:i5 Gvr for ncw Gold sample 
and SNLS data, respectively. These values are in agree- 
ment with the age of old stars [53] . Here we take three 
OHRG for comparison with the power-law dark energy 
model, namely the LBDS 53W091, a 3.5-Gyr old radio 
galaxy at z = 1.55 [56], the LBDS 53W069 a 4.0-Gyr 
old radio galaxy at z = 1.43 [57] and a quasar, APM 
08279 + 5255 at z = 3.91 with an age of t = 2.ll8;?Gyr 
[58]. The latter has once again led to the "age crisis". 
An interesting point about this quasar is that it cannot 



TABLE II. The best values for the parameters of 
f{R) modified gravity with the corresponding age for the 
universe from fitting with SNIa from new Gold sam- 
ple and SNLS data, SNIa+CMB, SNIa+CMB+SDSS and 
SNIa+CMB+SDSS+LSS experiments at one and two a con- 
fidence level. The value of Q m was determined according to 
equation (19) 



TABLE III. The value of r for three high redshift objects, 
using the parameters of the model derived from fitting with 
the observations at one and two a level of confidences. 



Observation 


X 




Age (Gyr) 


OIN Adl 11c W VjrUlLl 1 


6 2 07+ - 230 

°- zu ' -0.147 

6 207+ 445 

°- zu ' -0.515 


n 276+ - 376 

u - z ' D -0.240 

0.276±°;™ 


14 71+ ' 41 
ri- ' r_o.23 

i4.n±8:» 


SNIa(new Gold) + 
CMB 


6 190 +a224 

O.X»U_ 242 

fi iqn +0 - 433 
o.iyu_ 503 


0.278«; 36 7 ^ 
0.2781°,;™ 


14.69i°; 39 7 
U.69±° f 2 


SNIa(new Gold)+ 
CMB+SDSS 


fi 1fiq+ 0170 

°- iDy -0.184 

6-169™ 


0.28lt°;gl 
0.281+o-j 3 ? 


i4.56i8:» 

14.56lS;S 


SNIa(new Gold) + 
CMB+SDSS+LSS 


6.192l° ; 1 1 77 

fi 1 qo+0.321 
o.i»z_ 368 


n 070+0.273 

U.^( O-0.278 

n 070+0.524 

u - z '°-0.278 


14-691S:" 

i4.69l?:SS 


SNIa (SNLS) 


6.484i° ;»^ 
6.484i° ;^ 


0.2461°,;^ 

246+ 324 
u.zto_ 246 


13.69lo;l? 
13-691"; 37 


SNIa(SNLS)+ 
CMB 


6.477iH^ 

+0.194 
D -^' '-0.200 


246+ 162 

u - z ^°-0.165 
U.ZW_ 246 


13.681°;!? 
13.68l ;|l 


SNIa(SNLS)+ 
CMB+SDSS 


6.4281°;™° 
6-4281°;^ 


n oc;o+0.148 

u.zoz_ 152 
0.252«;^ 2 


13.60lg-.JS 
1 3 fifl +0,32 


SNIa(SNLS)+ 
CMB+SDSS+LSS 


fi 4 oo+0.089 

o.too_ 091 

64 oo+0.176 

o.too_ 183 


0.252t° : 1 1 47 
0.252_ 252 


13.6ll° ;i« 

i3.6ilS:S 



Observation 


LBDS 
53W069 
z = 1.43 


LBDS 
53W091 
z = 1.55 


APM 
08279 + 5255 
z = 3.91 


SNIa (new Gold) 


i oo+0. 05 
r-^J-0.03 

1 90+0.10 

1 - zo -0.09 


1 oo+0. 05 
r-J^-0.03 

i o 9 +0.11 
■■■•"■^-O-IO 


0.86lg; 37 
0-86lg; 3 8 8 


SNIa(new Gold)+CMB 


i 90+O.O5 
1 - zo -0.04 

i 90+O.IO 
1 - zo -0.08 


i oo+0. 05 
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1 09+O.II 
l-J^-0.09 


0.86lg; 37 
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SNIa(new Gold)+CMB 
+SDSS 


i 90+O.O3 
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i qo+0. 04 
r-J^-0.04 
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1 -' ;>z -0.07 
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0.85±g; 3 07 


SNIa(new Gold)+CMB 
+SDSS+LSS 


i 90+O.O4 
1 - zo -0.03 

i 90+O.O7 
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1 q 9 + 0.04 

1 - JZ -0.04 
i q 9 +0.08 
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SNIa (SNLS) 


1 1fi+ 002 
1 1fi+ 002 


1 91+O.O2 
1,z± -0.02 

i 91+O.O2 
1 " £ - L -0.02 
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be accommodated in the ACDM model [59] . To quantify 
the age-consistency test we introduce the expression r as: 

t(z;X ) t(z- 7 X )H 
t=— = — — s— , (60) 

tobs lots -"0 

where t(z) is the age of universe, obtained from the equa- 
tion (59) and t a b s is an estimation for the age of old cos- 
mological object. In order to have a compatible age for 
the universe we should have r > 1 . Table III reports the 
value of r for three mentioned OHRG with various ob- 
servations. We see that f(R) modified gravity with the 
parameters from the combined observations, provides a 
compatible age for the universe, compared to the age 
of old objects, while the SNLS data result in a shorter 
age for the universe. Once again, APM 08279 + 5255 at 
z = 3.91 has a longer age than the universe but gives 
better results than most cosmological models [60,61]. 

VII. CONCLUSION 

Here in this work wc obtained the dynamics of uni- 
verse with f(R) = \J R 2 — Rq gravity and compared it 
with recent cosmological observations. This compari- 
son has been performed with Supernova Type la Gold 
sample and SNLS supernova data, CMB shift param- 
eter, location of baryonic acoustic oscillation peak ob- 
served by SDSS and large scale structure formation data 
by 2dFGRS. The best parameters obtained from fit- 
ting with the new Gold sample data are: h = 0.63, 
X Q = 6.192t° Q {j 7 7 provides tt m = 0.278±g;|^| at lcr con- 
fidence level with x 2 min / N d.o.f = 0.900. Using the SNLS 
supernova data, the best fit value for model parameter 
is: X = 6.433±g;g|? gives O m = 0.252±g;^ at la con- 
fidence level with Xmin/Nd.o.f = 11.486. We also per- 
formed the age test, comparing the age of old stars and 
old high rcdshift galaxies with the age derived from this 
model. From the best fit parameters of the model us- 
ing new Gold sample and SNLS, we obtained an age of 
14.69±g;l| Gyr and 13.61±g;^| Gyr, respectively, for the 
universe which is in agreement with the age of old stars. 
We also chose two high rcdshift radio galaxies at z = 1.55 
and z — 1.43 with a quasar at z = 3.91. The ages of the 
two first objects were consistent with the age of universe, 
means that they were younger than the universe while 
the third one was not. 
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